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Temporal l/f a Fluctuations from Fractal Magnetic Fields 
in Black Hole Accretion Flow 
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Abstract 

Rapid fluctuation with a frequency dependence of 1//" (with a ~ 1 — 2) is characteristic of radiation 
from black-hole objects. Its origin remains poorly understood. We examine the three-dimensional mag- 
netohydrodynamical (MHD) simulation data, finding that a magnetized accretion disk exhibits both 1// Q 
fluctuation (with a~2) and a fractal magnetic structure (with the fractal dimension of D^l.9). The frac- 
tal field configuration leads reconnection events with a variety of released energy and of duration, thereby 
producing l/f a fluctuations. 
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Introduction 

Apparently random temporal fluctuations from Galac- 
tic black-hole candidates (BHCs; van der Klis 1995) and 
from active galactic nuclei (AGNs; Ulrich, Maraschi & 
Urry 1997) have led many astronomers to recognizing 
how complex the nature behaves. The light curves are 
neither periodic nor random around some mean. Rather, 
they are seemingly composed of shot events with a va- 
riety of peak intensities and durations (Negoro et al. 
1995). Number of analyses of X-ray light curves and op- 
tical AGN light curves show that power spectral density 
(PSD) is flat at lower frequencies (/) and is power-law 
(oc f~ a with a ~ 1 — 2) at higher frequencies. The break 
frequency corresponds to a few second for BHCs and to 
a few years for AGNs. More, l// a fluctuations are ubiq- 
uitous in natural behavior, although their origins have 
been unsolved (Tajima & Shibata 1997; Cable & Tajima 
1997). The significance of the 1//" noise is that it con- 
tains a long-term memory (Press 1978). It has been a 
puzzle how 1// Q fluctuations can arise in black-hole ac- 
cretion flows (or disks) under realistic circumstances. 

Among number of suggestions for a possible mecha- 
nism of variability, the most promising one is magnetic 
flares (Wheeler 1977; Galeev, Rosner & Vaiana 1979). It 
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has been established through the X-ray observations by 
the Yohkoh satellite that the solar flares are triggered by 
magnetic reconnection (Shibata 1996). In fact, solar soft 
X-ray variation exhibits l/f a fluctuations (UeNo et al. 

1997) . Similarly, sporadic magnetic reconnection events 
which can occur in accretion flow may be responsible for 
the variability of black-hole objects, as well (Mineshige, 
Kusunose & Matsumoto 1995). 

To prove this conjecture, we examine the three- 
dimensional data of global, MHD disk calculations first 
made by Machida, Hayashi & Matsumoto (1999). They 
calculated how magnetic field evolves in a rotating disk 
initially threaded by toroidal (B v ) fields. Since no cool- 
ing is taken into account in computations, the simulated 
disk is advection-dominated (Kato, Fukue & Mineshige 

1998) , rather than radiation-dominated as in the stan- 
dard disk. Then, the system we analyze corresponds to 
BHCs in the hard (low) state, in which fluctuations are 
largely enhanced and whose spectra can well be repro- 
duced by advection-dominated flow (Narayan, McClin- 
tock & Yi 1996). Magnetic fields are amplified with time 
via a number of MHD instabilities together with differen- 
tial rotation. The maximum field strength is determined 
either by field dissipation by reconnection or field escape 
from accretion flow via Parker instability. As a result, 
the mean plasma (3, the ratio of gas pressure to magnetic 
pressure, finally reaches ~10 irrespective of initial values 
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of 0. Locally, however, even low-/? (< 1) regions appear; 
inhomogeneous structure arises consequentially [see also 
similar discussion by Abramowicz et al. (1992) but for 
non-magnetic cases]. 

In the present paper, we thus analyze temporal and 
spatial behavior of the accretion flow in order to clarify 
the origin of fluctuations in black- hole objects. Re- 
sults of the analysis are presented in the following section. 
The final section is devoted to discussions. 



2. Temporal and Spatial Analysis of Accretion 
Flow 

Figure 1 displays the light curves (A) and their PSDs 
(B) of the simulated disk obtained in the quasi-stationary 
state. Here, we assume that magnetic reconnection 
events contribute much to the energy output and thus 
calculated time variation of J rjj 2 dV (with 77 and j being 
the electric resistivity and electric current density) inte- 
grated over almost whole disk. High frequency sides of 
the PSDs show nearly power-law decline with an index of 
a ~ 2, in agreement with the observations. This is, in a 
sense, amazing that without any fine-tuning of parame- 
ters or special assumptions an MHD disk naturally gives 
rise to l/f a fluctuations over three orders of frequency 
ranges. On the other hand, low frequency sides flatten at 
frequencies smaller than a reciprocal of several rotation 
timescale at a reference radius, which is also consistent 
with the observations (van der Klis 1995; Ulrich et al. 
1997). 

The appearance of 1// Q fluctuation, or more specifi- 
cally the presence of a long-term time correlation implies 
a long-distance spatial correlation in the distribution of 
magnetic fields. It is thus tempting to examine the spa- 
tial magnetic-field distribution. We specially pick up the 
quantity, j / p, the ratio of absolute value of electric cur- 
rent density to matter density, since it is a good indicator 
regarding a trigger of reconnection (Parker 1994; Ugai 
1999). In fact, it is shown by MHD simulations that a 
fast reconnection, as is observed in solar flares, occurs 
when the electric resistivity becomes anomalously high 
in localized regions (Tajima & Shibata 1997). Such a 
local, anomalous resistivity can be achieved where elec- 
tron drift velocity (which is proportional to j/p) exceeds 
a critical value (Yokoyama & Shibata 1995). Hence, any 
regions with high j/p values are all good candidates for 
a next reconnection site. 

We plot in figure 2 a snapshot of the spatial j / p distri- 
bution on a horizontal plane slightly above the equatorial 
plane. The panel roughly demonstrates to what extent 
reconnected area expands, once reconnection is initiated 
somewhere. We notice that the distribution is quite in- 
homogeneous; patchy patterns are visible everywhere in 
figure 2. Importantly, there seems to be no typical size 
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Fig. 1.. (A) Typical light curves of the simulated MHD disk. 
Here, we assume that radiation is predominantly due to field 
dissipation by magnetic reconnection, thus plotting the tempo- 
ral variation of the rjj integrated over almost whole disk. Elec- 
tric resistivity r\ is assumed to be r}§ X MAX[j'/p — (j/p) C rlt , 0], 
since magnetic reconnection seems to happen where local j/p is 
larger than some critical value. (B) The power-spectral densi- 
ties of light curves above. Dotted line is a least square fit with 
a broken power-law function for (i/p)crit of 200 (open circles). 
It is of great importance to note that general behavior does not 
depend on the values of (j//p) C rit- The unit of time is rotation 
period at a reference radius where the center of initial torus is 
located. The rotation timescale used here will correspond to a 
few seconds for BHCs and a few years for AGNs for in realistic 
situation. 



of each patch. The presence of fractal structure is sus- 
pected. 

To confirm this idea, we made a fractal analysis for 
the three-dimensional MHD disk with mesh point num- 
bers of (N x ,N y ,N z ) = (100,100,25); namely we first mark 
the sites where j/p exceeds some critical value, (j / p) C rit, 
and name any assembly of the marked sites clusters. We 
then count the numbers of clusters according to the clus- 
ter sizes (i.e., volume) for the data sets at five different 
timesteps and average the counts. The resultant distri- 
bution is plotted in figure 3 with thick lines. Surprisingly, 
the cluster size is distributed in a power-law fashion over 
three orders of magnitudes of cluster sizes, from a few 
larger clusters to numerous smaller clusters. 

For comparison, we calculated a random distribution 
in the following way: we evaluate first how much fraction 
of the entire volume is covered with the marked sites in 




Fig. 2.. Color contour map of the (j / p) distribution on a horizontal 
plane slightly above the equatorial plane. Here, values of (j/p) 
where colors change are as follows; 70 from white to blue, 120 
to green, 200 to yellow, and 300 to red, respectively. 



Fig. 4.. Relation between the volume (s) and mean radius (r) of 
each cluster. The lines of s oc r, s oc r 2 and s « r 3 are also 
depicted. The least square fit shows s oc r D with D ~ 1.9. 
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Fig. 3.. Histograms of clusters with j/p > (i/p) cr it = 200 (by the 
thick solid lines) and those with j/p > 400 (by the thick dotted 
lines) as functions of the cluster size (volume). Power-law dis- 
tribution is realized over three orders of cluster sizes; roughly, 
D(s) oc s~ 2 . For comparison, we also plot the same but for 
a random distribution (by the thin dashed lines). Obviously, 
large clusters are missing. 

figure 2, finding about 5.5% for j/p > 200. We next put a 
random number between and 1 in each site of the three- 
dimensional box, and mark the sites where the number 
exceeds 0.945. We then repeat the same procedure done 
for the MHD disk and plot the resultant histogram in 
figure 3 with the thin dashed line. Clearly, there are no 
very big clusters in the random model. In other words, 
a long-distance j / p correlation is lost there. 



Magnetic fields in the disk have fractal structure. To 
find a fractal dimension, we plot a size (volume) of each 
cluster as a function of its mean radius in figure 4. Here, 



the cluster mean radius is defined as r = 




— i"cm 
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s 

with rcM = 7^] r i f° r each cluster with a size of s(> 1), 

i 

where r, is the coordinates of the z-th site belonging to 
the cluster. From the fitting, we find roughly s a r fl 
with D = 1.9. 

3. Discussions 

Then, we address two key questions: what relates the 
fractal structure to the temporal fluctuation in the sim- 
ulated MHD disk? How can such a fractal distribution 
arise? Turbulence, in general, is known to exhibit fractal 
behavior (Procaccia 1984) and this fact may be related 
to our present finding, but predicted fractal dimension is 
D ~ 2.6, differing from that of the present case. Alter- 
natively, we note the notion of self-organized criticality 
(SOC; Bak 1996; Jensenl998), one of the most attractive 
concepts developed in the study of complex systems. 

Bak, Tang & Wiesenfeld (1988) proposed a sand-pile 
model to describe a system exhibiting l/f a fluctuation. 
Suppose that we fall sand particles one after another on 
a table. Fallen sand particles will form a pile, onto which 
another sand particle will be added. When a slope of 
the pile in either direction exceeds a critical value, an 
avalanche occurs and sand particles will slide down in 
that direction. Then, the system spontaneously evolves 
to and stays at SOC. In our case, addition of a sand 
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particle corresponds to energy input to magnetic fields, 
while the critical slope corresponds to the critical j/p 
over which energy dissipation occurs via reconnection 
(Mineshige, Takeuchi & Nishimori 1994). 

For systems in a SOC state, long-distance spatial com- 
munication among different sites is naturally built up and 
it yields long-term time correlation. If each flare light 
curve is expressed by a time-symmetric profile with expo- 
nential grow and decay about t = 0, L(t) oc exp(— \t\/r) 
with r (< r max ) being constant, its PSD is P T {f) oc 
e 2 /(l + 47t 2 / 2 t 2 ) 2 . If the energy (e) of each flare is dis- 
tributed as N(e) oc e~ p and each flare duration is related 
to energy as, £oct d , the total PSD becomes 



P (f) = £/M/Mt)At 

T 

r max t 2D 

K ^ (1 + 4tt 2 / 2 t 2 ) 2 



Hf) « ( J 



(3-p)D 



F(f) (Ar-0). (1) 
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slowly varying function of / for / > l/(27rr max ). We 
may regard e oc s (volume of a clump) and r oc r (mean 
radius of a clump). Then, we find from the simulation 
that p ~ 2 and D ~ 2 (see figure 4). Hence equation 
(1) leads P(/) oc /~ 2 , in agreement with the numerical 
result (see figure 1). That is the reason why the frac- 
tal magnetic field produces l// a fluctuations (Takeuchi, 
Mineshige & Negoro 1995; Kawaguchi et al. 1998). 

One of the most conspicuous natures of the SOC is 
its ubiquity; namely, it is supposed to describe various 
non-equilibrium open systems, such as earthquakes, for- 
est fires, evolution of biological species, and traffic flows 
(Bak 1996). In astrophysical context, it is important to 
note that coronal magnetic fields in the Sun are suggested 
to be in a SOC state (Lu & Hamilton 1991; Vassiliadis 
et al. 1998) as well, thus exhibiting power-law occurrence 
rate of flares and 1//" fluctuations in solar flare curves 
(UeNo et al. 1997). 

Gamma-ray bursts (GRBs) also occasionally exhibit 
1// Q fluctuations (Beloborodov, Stern & Svensson 1998). 
In some models of GRBs which involve merger of two 
compact objects (white dwarf, neutron star, and black 
hole) a sort of an accretion disk is thought to be formed 
by debris of one component around the other (Meszaros 
1999). The situation could be similar. We thus expect 
frequent reconnection events with a smooth size (ampli- 
tude and duration) distribution to occur in GRBs, which 
will give rise to 1// Q fluctuations (Panaitescu, Spada & 
Meszaros 1999). Likewise, any other magnetic systems, 
regardless of system size, may show similar effects. 
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